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Low-energy polaritons in semiconductor microcavities are important for many processes such as, 
e.g., polariton condensation. Organic microcavities frequently feature both strong exciton-photon 
coupling and substantial scattering in the exciton subsystem. Low-energy polaritons possessing 
small or vanishing group velocities are especially susceptible to the effects of such scattering that 
can render them strongly localized. We compare the time evolution of low-energy wave packets in 
perfect microcavities and in a model Id cavity with diagonal disorder to illustrate this localization 
of polaritons and to draw attention to the need to explore its consequences for the kinetics and 
collective properties of polaritons. 

PACS numbers: 78.66.Qn, 78.40.Me, 71.36.+C 



I. INTRODUCTION 

Planar semiconductor microcavities have attracted 
much attention as they provide a method to enhance and 
control the interaction between the light and electronic 
excitations. When the microcavity mode (cavity photon) 
is resonant with the excitonic transition, two different 
regimes can be distinguished based on the competition 
between the processes of the exciton-photon coupling and 
damping (both photon damping and exciton dephasing) . 
The weak coupling regime corresponds to the damping 
prevailing over the light-matter interaction, and the lat- 
ter then simply modifies the radiative decay rate and the 
emission angular pattern of the cavity mode. In contrast, 
in the strong coupling regime the damping processes are 
weak in comparison with the exciton-photon interaction, 
and the true eigenstates of the system are mixed exciton- 
photon states, cavity polaritons. This particularly results 
in the appearance of the gap in the spectrum of the exci- 
tations whose magnitude is established by Rabi (split- 
ting) energy. In inorganic semiconductors the strong 
coupling regime has been investigated extensively, both 
experimentally and theoreticallyfi^iii^&i and the dy- 
namics of microcavity polaritons is now reasonably well 
understood^ These studies continue to expand because 
of the prospect of important applications such as the po- 
lariton laser related to the polariton condensation in the 
lowest energy state, 5 ^ ^- 2 . 13 ! 14 ! 15 ! 16 ! 17 

In another development, organic materials have been 
utilized in microcavities as optically active semiconduc- 
tors. In many organic materials excitons are known 
to be small-radius states, Frenkel excitons, which in- 
teract much stronger with photons than large-radius 
Wannier-Mott excitons in inorganic semiconductors. 
The cavity polaritons, therefore, may exhibit much 
larger Rabi splittings on the order of few hundreds 
of mcV; large splittings have in fact been observed 
experimentally, 18 - 19 i 20 i 21 i 22 - 23 i 24 i 25 i 26 i 27 i 28 At the same 
time, Frenkel excitons typically also feature substan- 
tially stronger interactions with phonons and disorder - 



electronic resonances in both disordered and crystalline 
organic systems are frequently found rather broad and 
dispersionless. It is thus likely that manifestations of 
exciton-polaritons in organic microcavities could be quite 
different from the corresponding inorganic counterparts. 

In this paper we are concerned with the nature and 
dynamics of the low energy exciton-polaritons in organic 
microcavities, the states of particular importance for the 
problem of polariton condensation. Our goal here is to il- 
lustrate some qualitative features of the dynamics in both 
perfect and disordered systems and, thereby, to draw at- 
tention to the need of more detailed experimental and 
theoretical studies to elucidate conditions for the polari- 
ton laser operation based on organic systems. 

The bare planar cavity photons are coherent wave ex- 
citations with a continuous spectrum and whose energy 
e(k) = efc depends on the magnitude k = |k| of the 2d 
wave vector k: 

efe = (A 2 + ? l Vfc 2 / e ) 1/2 , (i) 

where A is the cutoff energy for the lowest transverse 
quantization photon branch we restrict our attention to, 
c is the speed of light and e the appropriate dielectric 
constant. 

In the vicinity of the excitonic resonance, w s (e be- 
ing the exciton energy), strong exciton-photon coupling 
leads to the formation of new mixed states of exciton- 
polaritons whose two-branch (E±) energy spectrum fea- 
tures a gap as, e.g., illustrated in Fig.[TJ Especially inter- 
esting are systems with detuning A — e small in compar- 
ison with the Rabi splitting. In the absence of the exciton 
scattering processes, cavity polaritons are also clearly co- 
herent excitations that can be well characterized by wave 
vectors k and have energies E±(k). Frenkel exciton scat- 
tering (due to phonons and/or disorder) in organic sys- 
tems with weak intermolecular interactions results in the 
exciton localization: excitons propagate not as coherent 
wave packets but by hopping; in such a regime the wave 
vector k is no longer a "good" quantum number. 
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FIG. 1: The energy spectrum of the polaritonic eigenstates 
in a Id model microcavity described by the Hamiltonian (|10[) 
with N = 1500 molecular sites and N = 1500 photon modes. 
Parameters of the systems are as follows: average exciton 
energy e = 2 eV, cavity photon cutoff energy A = 1.9 eV, 
dielectric constant e = 3, the exciton-photon interaction en- 
ergy 7 = 0.15 eV. The dashed lines show the energy eigenval- 
ues in the system without exciton disorder, a — 0, the solid 
lines correspond to the spectrum in the system with disorder, 
a — 0.03 eV. Here the energies are shown as a function of state 
"number" sorted in an increasing energy order, separately for 
the lower (LP) and upper (UP) polariton branches. In the 
system without disorder, the state numbers would be imme- 
diately convertible to the corresponding wave vectors. On this 
scale, the UP branches of perfect and disordered systems are 
hardly distinguishable. 



5k < k. 



(3) 



If 8E(k) is weakly fc-dependent, then, evidently from 
Eq. @, condition © would be necessarily violated in 
regions of the spectrum where the group velocity v g (k) 
vanishes. For the LP branch, as seen in Fig. [TJ this oc- 
curs both at its lower- and higher-energy ends. Reference 
|29| provided estimates of the corresponding "end-points" 
fc m i n and fc max , where 5k ~ k, for organic planar micro- 
cavities within the macroscopic electrodynamics descrip- 
tion of polaritons. It was then anticipated that exciton 
scattering would render eigenstates corresponding to par- 
ent polaritons with k < fc m j n and k > fc max spatially 
localized in accordance with uncertainty relations like 
8x5k x ~ 1. As we discussed above, at the higher-energy 
end (k > fc max ) of the LP branch, the eigenstates are 
practically bare exciton states in nature. At the lower- 
energy end (k < fc m ; n ), one would deal with localized po- 
laritons having comparable exciton and photon contents, 
particularly for the detuning |A — e\ small with respect to 
the Rabi splitting . Further numerical calculations^!^ 
for Id microcavities with diagonal exciton disorder con- 
firmed this qualitative picture for the polariton states. 
Below, we will use a similar Id microcavity model to il- 
lustrate the nature of the low-energy LP states as well as 
time evolution of low-energy wave packets. 



As we are actually interested in lowest energy polari- 
ton states, we will constrain further discussion mostly to 
the lower polariton (LP) branch. One can quickly notice 
that, since the photon energy ([1]) rapidly increases with k, 
higher-fc photon states would be only very weakly inter- 
acting with exciton states. Therefore a very large number 
of eigenstates of the system with energies close the upper 
end of the LP branch (see Fig.fl]) are essentially reflective 
of the bare localized exciton states with an incoherent 
propagation mode. For the remaining fewer and lower 
energy states of the LP branch, however, the exciton- 
photon interaction is strong and one is then faced with 
an interesting interplay of the bare localized nature of 
the "exciton part" of the polariton and the bare coher- 
ent character of its "photon part" . 

Transparent physical arguments were used in Ref. [2^ 
involving the indeterminacy (broadening) of the polariton 
wave vector k owing to the exciton dephasing from the 
general relation^ 



5E{k) 

hv g {ky 



(2) 



where v g (k) is the group velocity of parent polaritons in 
the perfect systen and 5E(k) is the energy broadening 
due to scattering. Based on Eq. @, one can at least 
distinguish polaritonic states where k is relatively well- 



II. DYNAMICS OF LOW-ENERGY WAVE 
PACKETS IN PERFECT MICROCAVITIES 

Before proceeding with a model analysis of polaritons 
in a disordered system, we will briefly discuss the time 
evolution of wave packets in perfect microcavities where 
all polaritons are coherent states well characterized by 
their wave vectors. Not only will this establish a compar- 
ative benchmark but is useful in itself as such dynamics 
reflects features of the polariton spectrum, and hence of 
the exciton-photon hybridization. 

Of course, specific features of the low-energy wave 
packets stem from the fact that the polariton dispersion 
near the bottom of the LP branch (k ~ 0) is manifestly 
parabolic: 



w(k) ~ lu + ak 2 , a = H/2M, 



(4) 



M being the cavity polariton effective mass, which makes 
the broadening of wave packets a relevant factor. Con- 
sider now a wave packet formed with the states close to 
the branch bottom: 



U(r,t) 



-iujot 



J dkA(k) e l ( k r - Qfc2t ). 



(5) 



It is convenient to choose the weight amplitude function 
A(k) Gaussian: A(k) = (/3/2tt 3 ) 1/2 exp [-(3(k - k ) 2 ] , 
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centered at wave vector ko. With this amplitude func- 
tion, Eq. (fSJ) yields 



|[/(r,t)| 2 = C(i)exp 



ft(r - 2ak t) 2 
' 2(ft 2 + a 2 t 2 ) 



(6) 



for the time evolution of the spatial "intensity" of the 
wave packet, 



C{t) 



ft 



2n(ft 2 + a 2 t 2 ) ' 



Equation ^ describes a Gaussian-shaped wave packet in 
2d whose center 

r c (t) = v g t, v s = 2ak , 

moves with the group velocity consistent with the disper- 
sion and whose linear width increases with time in 
accordance with the Id variance 



s(t) = (ft + a 2 t 2 /ft) 



1/2 



(7) 



The total energy in the wave packet is conserved: with 
our choice of the amplitude function, 



dv\U{v.t)Y 



1. 



Of course, Eq. ([6]) can be derived as a product of two 
independent Id normalized evolutions such as 



\U(x,t)f = y/C(F) 



exp 



j3(x - 2ak t) 2 
' 2(ft 2 + a 2 t 2 ) 



(8) 



which we will be relevant in our discussion of Id micro- 
cavities, in this case the Id packet amplitude function 



A(k) = (/3/27r 3 ) 1/4 exp [-ft(k - k ) 2 



(9) 



The spatial broadening ([7]) features an initial value of 
ft 1 / 2 and a characteristic time % = ft /a such that, at 
times t 3> tb, the variance grows linearly with veloc- 
ity Vb = a/3 -1 / 2 . To appreciate the scales, some rough 
estimates can be made. So for the effective polariton 
mass M on the order of 10 -5 mo (mo being the vacuum 
electron mass), parameter a = h/2M ~ 5 • 10 4 cm 2 /s. 
Estimates in Ref. [29| made with the Rabi splitting and 
detuning ~ 100 meV yielded for microcavities with dis- 
ordered organics fc min ~ 10 4 cm -1 . Then for the wave 
packets satisfying 1 < /3 1 / 2 fc m ij 1 < 10, the characteristic 
time would be 0.2 < tb < 20 ps and the corresponding 
velocity 5 ■ 10 8 > Vb > 5 • 10 7 cm/s. In our Id numeri- 
cal example below we will use the value of parameter ft 
within the segment just discussed. 

We note that by changing physical parameters of the 
microcavity and organic material, as well as conditions 
for the polariton excitation, one can influence the dynam- 
ics described above. One should also be aware that the 
evolution times are limited by the actual life times r of 
small wave-vector cavity polaritons. Long life times r on 
the order of 10 ps can be achieved only in microcavities 
with high quality factors Q = lot. 



III. TIME EVOLUTION IN A Id 
MICROCAVITY WITH DIAGONAL DISORDER 

Finding polariton states in disordered planar microcav- 
ities microscopically is a difficult task which we are not 
attempting in this paper. As a first excursion into the 
study of disorder effects on polariton dynamics, here we 
will follow Ref. [3l] to explore the dynamics in a simpler 
microscopic model of a Id microcavity. Such microcav- 
ities are interesting in themselves and can have experi- 
mental realizations; from the results known in the theory 
of disordered systems^ one can also anticipate that cer- 
tain qualitative features may be common for Id and 2d 
systems. 

The microscopic model we study is set up in the fol- 
lowing Hamiltonian: 

H = ^{e + e n )ala n +^2e k blb k 

n k 

+ 7 E Jjff { elkna + e" lfcna b{a n ) . (10) 

nk ' 

It consists of a lattice of N "molecular sites" spaced by 
distance a and comprises the exciton part (a n is the exci- 
ton annihilation operator on the site n) , photon part (bk 
is the photon annihilation operator with the wave vec- 
tor k and a given polarization) as well as the ordinary 
exciton-photon interaction. The cavity photon energy e k 
is defined by Eq. (jTJ) , e represents the average exciton 
energy, while s n the on-site exciton energy fluctuations. 
We will use uncorrelated normally distributed e n with 
the zero mean and the variance a: 



a 2 5 n 



(11) 



The exciton-photon interaction is written in such a form 
that 27 yields the Rabi splitting energy in the perfect 
system. We chose to use the same number iV of photon 
modes, the wave vectors k are discrete with 2ir/Na in- 
crements. Our approach is to straightforwardly find the 
normalized polariton eigenstates (i is the state in- 
dex) of the Hamiltonian (flTJ|) and then use them in the 
site-coordinate representation: 



*(n) = (* p (n), * e (n)) 



|*(n)| 2 = l, 



(12) 



where and ^> c respectively describe the photon and 
exciton parts of the polariton wave function and n de- 
notes the nth site. 

We have tried various numerical parameters in the 
model Hamiltonian with the results being qualitatively 
consistent; the parameters we exploit in this paper have 
been chosen, on one hand, to be reasonably comparable 
with the experimental data in the output and, on the 
other hand, to better illustrate our point within a prac- 
tical computational effort. It should be kept in mind 
though that we consider a model system and the numer- 
ical values of results may differ, likely within an order of 
magnitude, for various systems. 
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FIG. 2: Examples of the spatial structure of the the photon 
part |^ P | 2 of the polariton eigenstates in a Id microcavity 
with disorder, (a) Four states, shown by different lines, from 
the very bottom of the LP polariton branch with energies 
within the range of 1.76-1.77 eV (see the spectrum in Fig.[l}. 
The inset shows one of these states in more detail. Dots in the 
inset correspond to the sites of the underlying lattice, (b) One 
state with a higher energy close to 1.82 eV. The inset shows 
part of the spatial structure of this state in more detail. Dots 
correspond to the lattice sites and spatial oscillations of the 
wave function are clearly seen. 

The numerical parameters are indicated in the caption 
to Fig. Q] and have been used to calculate the eigenstates 
of the Hamiltonian (TTOj) with N — 1500 for a cavity of the 
physical length L = Na =150 /im and a small negative 
detuning (A — e) = —0.1 eV. Figure [1] compares the en- 
ergy spectra in the perfect microcavity and in the cavity 
with one realization of the excitonic disorder, Eq. (jlip . 
a = 0.03 eV. It is apparent that the effect of this amount 
of disorder on the polariton energy spectrum per se is rel- 
atively small, except in the higher-energy region of the 
LP branch where eigenstates, as we discussed, are prac- 
tically of a pure exciton nature. 

The lower-energy part of the LP branch, however, cor- 
responds to the polariton states *F (fT!?]) in which the 
exciton and photon are strongly coupled (7 = 0.15 > 
A — e| =0.1 eV) with comparable weight contributions 
in 'Fp and A dramatic effect of the disorder is in 
the strongly localized character of the polaritonic eigen- 
states near the bottom of the LP branch, as illustrated in 
Fig. [Ha) (needless to say that the same behavior is ob- 
served for the states near the bottom of the UP branch 
which we we are not concerned with in this paper). Of 
course, both the photon *F p and exciton ^ e parts of a lo- 



calized polariton state are localized on the same spatial 
scale, however the exciton part of the spatial wave func- 
tion is more "wiggly" reflecting more the individual site 
energy fluctuations. For better clarity, in both Figs. [5] 
and [3] we show only the smoother behaving photon parts 
'Fp. Panel (a) of Fig. [5] displays examples of I'Fpl 2 for 
four states from the very bottom of the LP branch in 
a realization of the disordered system that are localized 
at different locations of the 150 fim sample. The inset 
to this panel shows the spatial structure of one of these 
states in more detail; it demonstrates both the spatial 
scale I of localization in this energy range (I ~ 1 fjxa with 
the used parameters) as well as a "macroscopic" size of 
the localized state in comparison with the lattice spacing: 
I > a = 100 nm. 

The states at the bottom of the LP branch can be 
characterized as strongly localized in the sense of kl < 1 
where A; is a typical wave vector of the parent polariton 
states in the perfect system. This feature may be con- 
trasted to the behavior at somewhat higher energies and 
at higher k > fc m ; n of parent states, where the disorder- 
induced indeterminacy of the fc-vcctor becomes small sat- 
isfying Eq. ([3]) so that k would appear as a good quantum 
number. As is known^S however, the multiple scattering 
should still lead to spatial localization of the eigenstates, 
now on the spatial scale I such that H > 1. Panel (b) 
of Fig. [5] illustrates the spatial structure of such a state 
with much larger I than in panel (a). The inset to panel 
(b) shows that the wave function in this case, within the 
localization length, exhibits multiple oscillations with a 
period on the order of 1/fc, which produces the black 
appearance on the scale of whole panel (b) . 

Having all the eigenstates of the system calculated, 
we are now in a position to study the time evolution of 
an initial polariton excitation, which we choose in the 
form of a wave packet I'F ) built out of the low-energy 
polariton states \^} of the perfect system: 

|* ) = ^A i |^)=^B l |vI/ l ). (13) 

i i 

Polaritons in the perfect system are ordinary plane waves 
and we used a discretized analog of Eq. © for the ampli- 
tude function A,-, the result is a Gaussian-shaped wave 
packet as illustrated in Fig. [3] by the long-dashed lines 
for the photon part of the polariton wave function. Am- 
plitudes Bi in Eq. (j!3[) are, on the other hand, expansion 
coefficients of the same initial excitation over the eigen- 
states l^j) of the system with disorder. The time evolu- 
tion of the initial excitation in the perfect system is then 
given by 

|iF°(i)) = Y, A i e- lE * t/h I*?), (14) 

i 

while the evolution in the disordered system by 

\y(t)) = Y / B l e- lE ' t / n \y i ), (15) 
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FIG. 3: Examples of the the time evolution of spatially iden- 
tical wave packets built out of the polariton eigenstates of 
a perfect Id microcavity as in Eq. (|9]) with the parameter 
/3 1//2 = 5 /im. For panels (a) and (b), the initial packet has 
zero total momentum, ko — 0; for panel (c) the initial packet 
has a finite momentum determined by fco = 10 4 cm -1 . Only 
the photon part l^p] 2 of the polariton wave function is dis- 
played. The initial packets are shown by long-dashed lines, 
results of the evolution after indicated times t axe shown by 
solid lines for the disordered system and by short-dashed lines 
for the perfect microcavity (except in panel (a) , where the lat- 
ter practically coincides with the initial packet.) 



where E® and E\ are the respective eigenstate energies. 

Of course, the evolution of the low-energy wave packet 
(I14p in the perfect microcavity takes place in accordance 
with our continuum generic description in Eq. (jSJ) (bar- 
ring small differences that may be caused by deviations 
from the purely parabolic spectrum) . This is clearly seen 
in panels (b) and (c) of Fig. © where the photon part 
of |vE r °(i)) a t indicated times t is displayed by the short- 
dash lines: mere broadening of the wave packet with no 
momentum (fco = 0) in panel (b), and both broadening 
and translational displacement (ko ^ 0) in panel (c). On 
the time scale of panel (a), the \^°(t)) state has not prac- 



tically evolved yet from l^ ) and is not shown on that 
panel. 

The time evolution of exactly the same initial polari- 
ton packets is drastically different in the disordered sys- 
tem, the corresponding spatial patterns of the photon 
part of arc shown in Fig. [3] with solid lines. First 

of all, the initial packet is quickly (faster than a frac- 
tion of ps) transformed into a lumpy structure reflecting 
the multitude of the localized polariton states within the 
spatial region of the initial excitation. Note that in our 
illustration here we intentionally chose the initial am- 
plitude function ([9]) with the parameter /3 1 / 2 = 5 /im 
large enough for the spatial size of the initial excitation 
to be much larger than the size of the individual local- 
ized polaritons at these energies (compare to Fig. [5Ja)). 
Importantly, however, that, while displaying some inter- 
nal dynamics (likely resulting from the overlap of various 
localized states), this lumpy structure does not propa- 
gate well beyond the initial excitation region over longer 
times. This is especially evident in comparison, when the 
broadening and motion of the packets in the perfect sys- 
tem is apparent (panels (b) and (c) of Fig. [3]). We have 
run simulations over extended periods of time (~ 100 ps) 
with the result that l'F(i)} remains essentially localized 
within the same spatial region. Of course, some details 
of ^(t)) depend on the initial excitation - see, e.g., a 
somewhat broader localization region in Fig. [3Jc) for the 
initial excitation with an initial momentum correspond- 
ing to fco = 10 4 cm -1 but the long-term localization in 
the disordered system appears robust in all our simula- 
tions. It would be interesting to extend the dynamical 
studies with participation of higher energy states such as 
in Fig. [DJb), this is, however, beyond the scope of the 
present paper. 



IV. CONCLUDING REMARKS 

The nature and dynamics of low-energy cavity polari- 
ton states are important for various physical processes 
in microcavities, particularly for the problem of conden- 
sation of polaritons into the lowest energy state(s). As 
was demonstrated in Ref. l29l low-energy polaritons in 
organic microcavities should be especially susceptible to 
effects of scattering/disorder in the exciton subsystem. 
The problem of disorder effects on polaritons in organic 
microcavities appears quite interesting as organic mate- 
rials would typically feature both strong exciton-photon 
coupling and substantial static and/or dynamic exciton 
scatter ing . In this paper we have continued a line of study 
in Ref. |3l| to look in some more detail at disorder effects 
on polaritons in a Id model microcavity. Our numeri- 
cal analysis has brought further evidence that low-energy 
polariton states in organic microcavities can be strongly 
localized in the sense of I < A, where / is the spatial size of 
localized states and A the wave length of parent polariton 
waves. (We have also found indications of weaker local- 
ization at higher polariton energies in the sense of / 3> A.) 
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Our illustrations have included demonstrations of local- 
ization not only via the spatial appearance of polariton 
eigenstates but also via the time evolution of different 
low-energy wave packets. 

On the physical grounds^ one should expect that 
low-energy polaritons in 2d organic microcavities would 
also be rendered strongly localized by disorder, as it 
would also follow from the general ideas of the theory 
of localization. 33 Further work on microscopic models of 
2d polariton systems is required to quantify their local- 
ization regimes. 

The strongly localized nature of low-energy polariton 
states should affect many processes such as light scat- 
tering and nonlinear phenomena as well as temperature- 
induced diffusion of polaritons. Manifestations of the lo- 
calized polariton statistics (Frenkel excitons are paulions 
exhibiting properties intermediate between fermi and 
bose particles) in the problem of condensation also ap- 
pear interesting and important. 

We note that one can exercise an experimental con- 
trol over the degree of exciton-photon hybridization and 
disorder by modifying the size of microcavity for various 
organic materials making such systems a fertile ground 
for detailed experimental and theoretical research into 



their physics. 

And the last remark. While we have specifically dis- 
cussed exciton-photon polaritons in organic systems, it 
is clear that some aspects have a generic character and 
could be applicable to other systems. This, for in- 
stance, concerns inorganic semiconductor microcavities. 
Both exciton-photon coupling and magnitudes of disor- 
der there, however, are much smaller, which would likely 
make any localization effects relevant only at very low 
temperatures. As one of recent examples of a very dif- 
ferent kind of systems, we will mention hybrid modes 
in chains of noncontacting noble metal nanoparticles 
where the interaction of photons and nanoparticles lead 
to plasmon-polaritons. 34 
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